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I. INTRODUCTION 

Firehose (FH) and Mirror Mode (MM) instabilities are 
both pressure anisotropic instabilities, which mainly hap- 
pen in high beta plasma, and have many applications 
in space and astrophysical physics. One (FH) happens 
when parallel pressure exceeding perpendicular one, and 
the other (MM) happens when perpendicular pressure 
exceeding parallel one. Both of these two modes have 
been discussed in literatures by many authors. The ba- 
sic properties of them can also be found in textbooks or 
monographs, see e.g., [8] and [30] . 

Mirror mode (MM) instability was identified and dis- 
cussed since [53], where the MHD theory is used, which 
is only suitable for long wavelength limit. [29] gives a 
kinetic description, and shows that this mode is not a 
simply fluid instability. [6] put forward to both the ef- 
fects of finite Larmor radius (FLR) and nonuniform, i.e., 
drift-mirror mode. j27j gives a discussion of the physical 
mechanism of the linear mirror instability in cold electron 
temperature limit. Recently, series papers by Pokhotelov 
et al addressed many details of MM via both fluid theory 
and kinetic theory, such as finite electron temperature 
effects [14] , non-Maxwellian velocity distribution [16] , fi- 
nite ion-Larmor radius wavelengths [TS] [H] . Gyrokinetic 
theory and simulation are introduced to discuss linear 
and nonhnear MM by Qu et al, [Mf^ . 

Firehose instability is another well-known plasma in- 
stability, which is in fact from the same branch of the 
kinetic Alfven wave (KAW), hence, we will not distin- 
guish them in the following sections. When isotropic but 
with FLR, we get KAW [7j^; while, when anisotropic 
but using long wavelength limit, we get the classical FH 
[31 [131 [21]. It is found by [31 (also a review paper [32]) 
that a significant new effect which had been neglected in 
the past will be brought in when allowing the ion gyro- 
radius to be finite. [3] extends [31] to discuss electron 
temperature anisotropy effects. 



FH and MM have been discussed in almost the same 
approaches due to their similarities, e.g., CGL fluid, ki- 
netic corrections. However, most the analytical solutions 
are reduced to just including one of them, and then, there 
are few discussions of coupling effects. For examples, 
when mirror mode is unstable, where is the FH (KAW) 
solution? Can FH be also excited to unstable? Can FH 
bring MM from a pure damping/growth mode to with 
real frequency? [26 just discusses the nonlinear growth 
of FH and mirror fluctuations, mainly via simulation. [5] 
has discussed the coupling indeed, but uses the fluid the- 
ory then not general. The gyrokinetic treatment of the 
stability of coupled Alfven and drift mirror modes in non- 
uniform plasma can be found in [111 I12j. however, where 
FLR effects are neglected. 

Here, we talk about linear coupling effects to answer 
the above questions. Firstly, a general anisotropic (bi- 
Maxwellian) 3-by-3 dispersion relation matrix is derived 
in gyrokinetic framework (then, with FLR effects but 
only for k\\/k±^ <^ 1) in Secjilj In Sec ]III[ we show the gy- 
rokinetic dispersion relation can reproduce the FH and 
MM solutions. In Sec IV using only cold electron as- 
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sumption Te/T^ <C 1, the matrix is reduced to 2-by-2 
for discussion the coupling of FH and MM. The dis- 
persion relation is solved both analytical and numer- 
ical, and is also compared with the full-kinetic code 
WHAMP (Waves in Homogeneous Anisotropic Multi- 
component Magnetized Plasma, [13]) • General stability 
properties for arbitrary Larmor radius at the (/3||,/3_l) 
plane are made clearly from both analytical and numeri- 
cal solutions and are also confirmed by Nyquist analysis 
in Sec]V| A summary is drawn in the last section. Ap- 
pendix gives some detailed discussions of the gyrokinetic 
dispersion matrix. 



II. GYROKINETIC DISPERSION RELATION 
MATRIX 

Here, we use the gyrokinetic (GK) theory [T] instead 
of the full kinetic (FK) theory to calculate the general 
dispersion relation matrix, which excludes the high fre- 
quency (r^ci ) modes automatically. The linear GK or- 
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dering is 

1. Small amplitude: Sf/F^ - e5ip/T - 5B/Bq - (5 < 
1; 



and, 



2. Low frequency: uj/^li ^ S; 

3. Anisotropic; k\\/k±^ ^ Pi/L^ 



5, k±pi 1. 



For detailed discussion of the gyrokinetic assumptions, 
one can refer to [1] or [10]. We comment here, the same 
ordering S used here for all small variables is just for 
convenience, and the third assumption k±pi ~ 1 is not a 
must. 

Electrons and ions are assumed both bi-Maxwellian 
distribution, 



Fn = 



1 



cxp 



Q:_L,|| = 



2r, 



1/2 



(1) 



Using the framework of _2J, we write the linearized 
equations of quasi-neutrality condition, vorticity equa- 
tion (or parallel Ampere's law) and perpendicular Am- 
pere's law to matrix form 







Si! 








L Bo 9i J 





css CAS Cms 

CSA CAA CmA 
CSM CAM CmM 



C _ = \ CSA CAA CMA ^ =0. 



5i! 

Bo 9. 



(2) 

Here, Coulomb gauge V • A = is used. S<j) is the per- 
turbed electrostatic potential, Sip = SA\\ -oj/ck^^ is a quan- 
tity related to the parallel component of the perturbed 
magnetic vector potential, and 5B\\ is the parallel compo- 
nent of the perturbed magnetic field. The final dispersion 
matrix can be gotten as. 



Til 



7ii. 



T±' Tile 



Tj_i ^ SnnpT 

rp T P — n2 ' 



, Kp' -2_ 
— — , — 



klv\'' 



5(t)\\ =5(t>-5ip {E\\ = -ifc||(50||), 

Toib) = /o(6)e-^ 

ri(6) = [Io{b) - hib)] e-^^2(6) = 2T,{b). 

where /„ is the first kind modified Besscl function. The 
dispersion relation is 



det C =0 



(4) 



We see from ([3|, the matrix is symmetrical, i.e., cas = 
CSA, CAM = Cm A and cms = csm, which is similar with 
the full- kinetic one (see e.g., Plj). When isotopic (yy^ g = 
1), Q reduces to exactly the same result as [TDj. At 
appendix, we also show that when drop small terms from 
reduced full-kinetic result in [31], it also gives the same 
result as this gyrokinetic one. 

We can see from ([3| , the solutions for ui is independent 
of which means always cx fcy. This is the main 
drawback of Q because of that fcy /k± is taken as small 
term and dropped, then not suitable to discuss the fine 
structures of the dispersion relation, e.g., the maximum 
growth rate 7max of mirror mode. However, one can also 
add the small terms when necessary by comparing with 
the full-kinetic result (see appendix). For example, [20] 
adds a k\\/k^_ term in integral form when discuss mirror 
mode. 



css Cas Cms 

C = I CSA CAA CMA 

_ csM Cam cmm 

-Ai + A2 -??,;A5 -A3 - A4 

-7?jA5 ^(1 A) - TyiAs TyjAe 
—A3 — A4 TjiXQ -1^ — A7 — Ag 



(3) 



where, 



A=^(l-ry.)[l-ro(fe.)] 
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(l-r;,)[l-ro(&e)] 



Ai = [1 + ^^Z{i,)Tp{bi)] + T [1 + ee^(Ce)ro(6e)] 

A2 = (1 - m) [1 - ro(6.)] + r (1 - 77e) [1 - ro(6e)] 

A3 = (l-??.)ri(6,)-^(l-r7e)ri(6e) 

A4 = 6^(C^)ri(6.) - ^ee^(ee)ri(6e) 

A5 = [l-ro(6,)]+ri[l-ro(6e)] 

As = ri(6,)-ri(6e) 

A7 = (1 - m)T2{h) + i(l - 7?e)r2(6e) 
As = £,rZ{i,)V2{b^) + ^^e^(Ce)r2(6e), 



III. SOLUTIONS OF FIREHOSE AND MIRROR 
MODE 

In ([3]), when only considering one matrix element, 
Css — gives ions sound wave (ISW, only consider 5cf)\\), 
CAA — gives shear Alfven wave (SAW, only consider 
Sij)) and Cmm — gives mirror mode (MM, only consider 
5B\\). Cap{a, (3 = S,A, M) is the coupling of each waves. 

A. Firehose (KAW) Branch 

As mentioned above, caa — represents the firehose 
and KAW solution, 



= f (1 + A) = 



[i-ro(6.)] + ri[i-ro(6e)] 



1 + ^ (1 - v^) [1 - ro(fe,)] + [1 - ro(^'.)] \ -(5) 

When small Lamor radius 6, ^ ^ Ij ro(^) — 1 ~ ^ " 
36^/4 ~ 1 - 6, ([5| reduces to. 



Pi 



cD^ = l + A = H-^(l-^,) + ^(l-^,), (6) 
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which is the classical firehose solution, and when 1 + A < 
0, this mode is unstable. When isotropic rji^e = 1; ^ = 0, 
^ reduces to 



L0^ = 



[i-ro(6,)] + T^[i-ro(6e)]' 



(7) 



Still using the small Lamor radius assumption to go fur- 
ther, Q reduces to 

^ 1 + 36,/4, (8) 
^ is the classical KAW solution [312]. 

B. Mirror Mode Branch 



At this case, cmm — 0, gives 

At - As = 0, 



(9) 



And, dropping the contributions from electrons (cold 
electrons assumption) in (|9|, gives, 



(10) 



(10 1 is the mirror mode solution. [20] adds a ki\/k± term 
to(10), then can be used to discuss the 7max of mirror 



mode, which give the same result as [TH] from full-kinetic 
theory. 



IV. COUPLING OF FIREHOSE AND MIRROR 
MODE 



In the above section, we use only one matrix clement 
of the 3-by-3 general dispersion matrix ([3|, then get ki- 
netic version of both FH (KAW) and MM, which shows 
accurately enough to discuss single mode by compared 
with full-kinetic results from previous authors. 

At this section, we solve the 2-by-2 matrix to discuss 
coupling effects of FH and MM. The matrix ([3| is general 
for arbitrary 6^ ej A.ej Ci,e under the gyrokinetic order- 
ing. When cold electrons {T^/Ti <^ 1), we have be ^ 1, 
T,T± > 1, css > csA,csM, then (5(/)|| should be very 
small, which means the parallel electrical field is short cir- 
cuited (few percents of cold electrons in density is enough, 
see [5j). We use this only assumption to reduce the 3-by-3 
matrix to 2-by-2 for discussion FH and MM coupling. 

Q reduces to. 



2^(1 -f A) 



2ru_ 

0±, 



— A7 — Ag 



(11) 



where, 
A 



^(l-,7.)[l-ro(6.)] + %^(l-^e) 



Pi 



(1-^0(1 



A5 ~ 1 - To{b,) 



ri(6. 



1 



As = 2C,Zi^,)T,{b, 



A7~2(l-7/,)ri(6,), 

' ^''''-ee^(ee)ri(&e) 



2^,Z{^i)T^{bi). (12) 



The small bi expansion is also written in ( 12 ) but only 
used when necessary. 



has given a very similar 2-by-2 matrix as (11) 



via full-kinetic theory, which contains also the small 
fc||/fc_L correction. But, it solves only MM by drop- 
ping the firehose correction in the mirror approximation 
\uj\ <^ A fluid theory is used to discuss the in- 

stability of mirror mode and firehose in [5], and several 
stable and unstable regions are given in different param- 
eters space. Since there is no kinetic effect, the results 
of [5] will be unreasonable. Ignoring FLR effects, [Ij] 
gives the gyrokinetic treatment of the stability of cou- 
pled Alfven and drift mirror modes in ID non-uniform 
plasma. 

For analytical tractable, we solve ( [TTj ) in long wave- 
length approximation in this section. This means bi <^ 1. 
Keep 0(1) -I- 0(6), (In]) changes to 



f3i 



[(l-r;,) + e.^(eO] 



1 



(1 + A)-1 



l^-f [(l-r;,)+e.^(CO] + 
[{l-r,,)+^,Zm [-=;^(H-A)-1] 



6(13) 



FH to MM. 



We use ( 13 ) to see the correction from MM to FH and 



A. Mirror mode correction to Firehose (KAW) 

Look the small FLR correction to firehose solution 
[il + A)/cj^-l]^0{b,), 



(1 + A) 

^ {ljgT-|[(i-»K)+g.^te)]~f>7.}b. 



(14) 



1. Cold ions 

When cold ions, |w| » |fc||Uj|[|, £_i 
6^(6) ^ -1 - 1/(2^2) -1, 



'^/\/% > 1. 
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-,2 ^ (1+A) 

=1, ;9, i<l 



1 + |6, (KAW) 



Z. Hot ions 



(15) 





ID aj zi a 



n a a im 



D iDSXiDsainiaaiim 



Here, we try to see that, when classical FH (KAW) 
is stable ( 1 + A > 0), whether MM can excite it 
to unstable. When hot ions, |w| ^ = 

^V^Vl + Asgn(d;) / y^^, (Note: in fact, w is complex) 



(1+A) 



UJ 



[A- 



3^(A- 



a-is5n(cj)6 (a - 0(1), 5-0(6;)) 



where, 



A; 



v/ttVI + Asgn(w) 



(16) 



(17) 



Due to — sgn((I)), and provided, a > 0, the solutions 
are always damped. 

To find a unstable mode, which means a < 0, 



a < 0^(1-^)^2 - 
4 



>o 



^6,77,A+(1-:^)S2<0, (18) 
^ ^ ^ ^4 ^ 

<0 >0 



(181 can only be satisfied when A < 0, i.e., mirror mode 



is unstable. While, if there is a domain satisfy both a < 
and A < 0, there must be at least one point satisfy both 
a = and A = 0, 



a = 0, 
A = 0, 



(1 



A = 0, 



0, 



5 = 0, 
A = 0, 



A = 0, 
- (1 - '7^ 



0, 



(19) 



( |19| can never be satisfied, then, there is no unstable 
domain, providing we use small argument expansion of 
ro,i(6,) and Z(e.). 

We can also contour plot (see FIG[T]) the a = bound- 
ary and A = Q boundary to show whether there is a 
region satisfy both (181 and A <Q. However, these con- 



tour plots can also help to show the sizes of each kind of 
ranges. 

We see that, indeed, none of the parameters (/3||i, f5±i) 
satisfy both a < and A < 0, which means that the 
unstable domain, if it exists, should be very very narrow. 
The only possibility is when To.! (6,) or Z{(^i) cannot be 
expanded, i.e., &i 1 or — 1. In next section, we 
will use Nyquist technique to prove that there is indeed 
no unstable domain for exciting the classical stable FH 
(KAW) in arbitrary parameters. 



FIG. 1: Contour plot of a = and A^Q. 



3. Numerical solutions 



Numerical solutions of (111 at FH (KAW) branch are 
jiven below, 





(a) 




-1.3 
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li 



(b) 



FIG. 2: Numerical solution and analytical solution ( 16 1 
firehose (KAW). 



FIG[2] shows, the analytical and numerical solutions 
are agreed very well. The coupling of MM is weak, and 
mainly brings KAW (FH) a damping. 
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B. Firehose correction to mirror mode 



At this case 



— A7 — As 



r(l + A)-Aj 



(20) 



We define Wq be the traditional mirror mode solution 
(zeroth order), i.e., — cjo/-v//3||i ^ 1, and 



^io^(Cio) — 



%-(l-77,)ri(fe,)/?. 



/3±»ri(fe,) 



(21) 



We can rewrite \20\ as, 
C.^(6) 



- (1 - 



r?i[ri(6i)-l]' 



2{t(l + A)-[l 



which brings just a very small correction and MM is still 
pure growing or damped. This means the coupling from 
KAW (FH) is also weak, and the traditional treatment 
of MM reasonable. To go further, we use the expansion 
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i^fnuij A//3||i. The ana- 



lytical solution in F IGp3| is by solving (21) and (22) with 
6Z(6) ~ i^/TTUj / ^ f3\\i, and then adding the correction 

~ 



term 
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Numerical solutions of (111 at mirror 

mode branch are also given. 

From FIGj3];b), we can see, for mirror mode, the an- 
alytical solution is only suitable when 6j ^ 1, because, 
when b, >- 1, ii = 1 . At FIG|3];c), the mir- 

ror mode assumption is totally broken when hi is large. 
WHAMP result is also shown. 



C. Physics Mechanism 

1. Mirror mode to firehose (KAW) 
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We find from the analytical and numerical solution 
above, that, the correction from mirror mode to firehose 
(KAW) is always bring not a unstable but a damping. 
This seems mean that the mirror cannot excite firehose 
but, on the contrary, absorb energy from firehose. 

It's widely accepted that the mechanism of mirror 
mode unstable is both the Landau type (wave-particle 
resonant) and the anisotropic free energy type. Then, it 
is physical reasonable that when mirror mode is unstable, 
it absorbs energy from firehose and cause which damped. 



(c) 



FIG. 3: Numerical solution of and analytical solution 
from (221, mirror mode. 



the free energy of non-Maxwellian distribution cannot 
compete to Landau damping). 



2. Firehose to mirror mode 

Mirror mode is always pure growing or damping 
(Re(|d;|) <C Im(|w|)), even though with firehose coupling. 
There is a narrow range (/3j_i >~ /3||i)i both mirror mode 
and firehose are damping, which may be caused by Lan- 
dau damping, similar with beam-plasma system (a range 



V. NYQUIST STABILITIES ANALYSIS 



The dispersion relation (11) can be rewritten as 



D{uj) = {Cj' -ao){ai-i,Z{U))+ uj'c' = 0, (23) 
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where 



(a) Nyquist Diagram oiD(m), ^.^^=2 5, p.j_=4, b.=D 25 



(b). Zoom m the box region of [a] 



1-Toib,) 



(1 + A),ai 



/9±^^l(6,) 



(1 



2 ^ ^4^1(6,) ~1] 

2ri(&,)[i-ro(6.)] 



Classical results of firehose and mirror mode instabili- 
ties, 



Firehose : 



Mirrormode 



stable : ao > 0, 
unstable : ao < 0. 

stable : Qfi > 0, 
unstable : ai < 0. 



(25) 



For different (/Sy , P±), there are three regions: (f ) ao > 
0, ai < 0, mirror mode unstable region; (2) < 0, ai > 
0, firehose unstable region; (3) ao > 0,ai > 0, both 
stable (damping) region. 



A. Mirror mode unstable region 

We use Nyquist technique to prove that when ao > 
0, ai < O(mirror mode unstable region), there is only 
one unstable mode (i.e., mirror mode). 



(26) 



D{u!) = (w^ - ao)(ai - ^iZir) +ui^c'^ - i(w^ - ao)^^Zu, 
When \uj\ oo, 

D{uj) ^ ui^{ai + l + c^) >0, 
Along — cxD < < oo, at ImD(cL') = 



(27) 
(28) 



— aQ,D — aQ(? > 0, 
(Ii 0, -D — — Ofgai > 0. 



The Nyquist diagram will draw as FIG|4j 



(29) 
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FIG. 5: Nyquist diagram of ( 23 1 verified by numerical, mirror 
mode unstable region. 



B. Firehose unstable region 



At this case ao < 0, ai > 0. When oo, we still 



have (28 1. And, along —oo < uj^ < oo, \v&D(uj) = at 



only w = 0, which gives D = — aoai > 0. Then, after 
drawing the Nyquist diagram, we can also get only one 
unstable mode (i.e., firehose). A result is shown in FlGj6] 



(s), Njquisl Diggram of D(ffl), P,y/=4.5, P,^=3, b.=0.25 




(b) Zoom in the boj region of [a] 
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Re(D) 
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FIG. 6: Nyquist diagram of ( |23[ ) verified by numerical, fire- 
hose unstable region. 





FIG. 4: Nyquist diagram of (231, mirror mode unstable 
region. 

We see from FIG|4j the mapping in complex D-plane 
indeed encircle the origin point, and one time, which 
means there is indeed only one unstable mode. A nu- 
merical result of the mapping (partly) is shown below in 
FlGEl 



C. Both stable (damping) region 

At this case ao > 0, ai > 0, the analysis is similar 
with the mirror mode unstable case, except that in Fig4, 
the point C > in .D-plane should be changed to C < 0. 
Then, the mapping in complex D-plane will never encircle 
the origin point, which means all the solutions are stable. 
A result is shown in FIG 17] 



VI. SUMMARY 

In this paper, we firstly give a general 3-by-3 gyroki- 
netic dispersion matrix for anisotropic bi-Maxwellain dis- 
tribution plasma to discuss low-frequency hydromagnetic 
waves, which is suitable for arbitrary &i^e, /?i,e, Cj,e under 
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[a], Nyquisl Diagram ^.^^=2 5, p.^=3, b.=0.25 




(b) Zoom in the bos region of (a) 




where, 



-0.5 D 0.5 1 15 2 
Re(D] 



FIG. 7: Nyquist diagram of ( 23 1 verified by numerical, stable 
region. 



the gyrokinetic ordering. This gyrokinetic matrix is use- 
ful to discuss various kinetic corrections to hydrodynam- 
ica waves. Some detailed solutions are given in [10 for 
isotropic plasma, e.g., analytical form of high /? KAW. 
However, the gyrokinetic matrix can also be extended to 
arbitrary distribution and arbitrary species, and also to 
nonuniform plasma. One can do this by following the 
framework in 15], see e.g., |12j for ID nonuniform plasma 
case. As an application example, we use this disper- 
sion matrix to discuss the coupling of firehose and mirror 
mode. At the cold electron assumption, the matrix is re- 
duced to 2-by-2, for the coupling of firehose (KAW) and 
mirror mode. The dispersion relation is solved both an- 
alytical and numerical, and shows consistent very well 
between each other. The results by the kinetic code 
WHAMP with all kinetic effects also show the reduced 
dispersion relation is appropriate. We find there exists at 
most one unstable solution at the {/3\\,/3±) plane for ar- 
bitrary Lamor radius, and the plane can be divided into 
three regions: firehose unstable region; mirror mode un- 
stable region; both modes stable (damping) region. The 
Nyquist analysis confirms these three regions. 
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VIII. APPENDIX 

At this appendix, we compare the full-kinetic disper- 
sion relation from [31] with our gyrokinetic one. 

Drop the electrons term as |31| by assuming be <ti 1, 
Q can be rewritten to, 

[t],DAC + v^V^E^A - ri,DB^]uj^ + ?7,6,(I + A) 

Cii^^D -A) + {rj,E - B f] = 0, (30) 



A=^(l-ry,)[l-ro(5,)]-|-^(I-,7,) 
^ = 2 ^ — 2 — 



c = 



2r;,ri(50 + 2^ + ^r2(5^ 

T_L 2 



Tie 2 



7?eTj_ 2 

[l-ro(5,)] - Ag, 



Dropping the 0{S^) term /c^w^/rif in keep only 
0(1) term in r\ and 77', the full-kinetic dispersion relation 
in 1311 can be written as 1331. 



-Ae\-) +(2A-M,)(1 + A)J^(-) 



+ -^(1 + A)^(1 + A) [(M.-^)] 



6,0(52) 



(1 + A)fe,7y, 



E - 



B 
2^ 



(32) 



where, 

, 1 - he-^' 



2A. E ^^ 2 ' - Oil) (33) 

n=l 



Dropping the high order 0{d ) terms in (32), we get 



exactly the same result as the gyrokinetic result ( 30 1 . We 
find in (32 1 when 6j ~ fcyU^/ri^, the biO{S'^) wiU jump 

to 0(i5^)!! So, we cannot drop the 0(5'') in (32) when 
bi^O. This is what [SI] found. 

In [31], the classical firehose solution is found in full- 
kinetic at the 0{6'^) terms, because biO{5^) terms can be 
dropped when 6i — J' ^ 0((5^). While in gyrokinetic, 
the solution is found in biO{5'^) terms, because we have 
ruled out the 0{6'^) terms by gyrokinetic assumptions. 
It should be just a coincidence that they give the same 
firehose result. The different can also be explained in 
another way: to get the classical firehose solution, the 
full- kinetic approach of [31^ uses the assumption bi — 
<C 0{5^), while the gyrokinetic approach uses b^ — > 

» 0(52). 

To discussion the fine structure of the hydrodynamica 
waves and instabilities, (32) can be seen as an extent 



version of our gyrokinetic dispersion relation, and won't 
go so complicated that can only be solved via numerical. 
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This is why we indeed have general dispersion relation 
for arbitrary distribution and with all kinetic effects in 
[28j and also can be numerical solved in very general form 
by WHAMP ("27, but should note that WHAMP is not 



suitable for strong damped waves, because the plasma 
dispersion function in that code is expanded via Pade 
approximation), but we still derive many other reduced 
dispersion relations to meet special desires. 
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